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Preparation, manipulation, and detection of strongly correlated states of quantum many body 
systems are among the most important goals and challenges of modern physics. Ultracold atoms offer 
an unprecedented playground for realization of these goals. Here we show how strongly correlated 
states of ultracold atoms can be detected in a quantum non-demolition scheme, that is, in the 
fundamentally least destructive way permitted by quantum mechanics. In our method, spatially 
resolved components of atomic spins couple to quantum polarization degrees of freedom of light. 
In this way quantum correlations of matter are faithfully mapped on those of light; the latter can 
then be efficiently measured using homodyne detection. We illustrate the power of such spatially 
resolved quantum noise limited polarization measurement by applying it to detect various standard 
and "exotic" types of antiferromagnetic order in lattice systems and by indicating the feasibility of 
detection of superfluid order in Fermi liquids. 



Introduction Future applications of quantum physics 
for quantum simulations, computation, communication, 
and metrology will require an extremely high degree of 
control of preparation, manipulation, and, last but not 
least, detection of strongly correlated states of quantum 
many body systems. Ultracold atoms offer an unprece- 
dented playground for realization of these goals. Sev- 
eral paradigm examples of strongly correlated states have 
been successfully realized, such as the Mott insulator, the 
Tonks gas, and the Bose glass (for a review cf. [1]). A 
standard way of analyzing such systems is by releasing 
the atoms from the trap and performing destructive ab- 
sorption spectroscopy, which only allows to measure the 
column density of the expanded cloud. Considerable at- 
tention has been thus devoted recently to novel meth- 
ods of detection, that allow for measuring (spin) density- 
density and other higher order correlation functions. One 
of those methods is atomic noise interferometry [2j , whose 
power is well illustrated in the recent observation of the 
bosonic and the fermionic Hanbury Brown-Twiss effect 
[3J 0] . Direct atom counting is another way to measure 
this effect, and to even go beyond it [5]; it can be realized 
directly with metastable Helium atoms [51 [7], or by using 
methods of cavity quantum electrodynamics (QED)|8J. 
Cavity QED is also essential in the recent proposals of 
Ref. [9j [10] , while Ref . [UJ proposes how to prepare and 
detect magnetic quantum phases using superlattices. All 
of the above approaches are, at least in some respects, 
destructive and frequently suffer from undesired atom 
number fluctuations inevitable in the preparation of the 
quantum states. 

Here we propose a unique method that allows for spa- 
tially resolved quantum non-demolition (QND) [12] de- 
tection of quantum states of ultracold atoms with inter- 
nal (pseudo-)spin degrees of freedom. Our approach is 
based on the idea of quantum noise limited polarization 
spectroscopy, demonstrated in Ref. [13 . Polarization 
spectroscopy itself consists in shining a polarized probe 



laser beam through the atomic spin system, and has been 
a subject of studies for years [2]. What is novel in quan- 
tum noise limited schemes is that quantum polarization 
degrees of freedom of light couple to the atomic spins, and 
in this way quantum fluctuations of the spins (magnetiza- 
tion) are faithfully mapped onto those of light; they can 
then be efficiently measured using homodyne detection of 
the transmitted probe. With an appropriate choice of pa- 
rameters, and provided that shot-noise-limited detection 
of light is achieved, this approach has proved very suc- 
cessful for QND-based quantum interfaces between light 
and atoms. In particular, atomic squeezing [15;, atomic 
entanglement, quantum memory, and teleportation have 
been achieved in this context (for a review see [16] and 
references therein). We have recently proposed to apply 
this approach to detect magnetic order of weakly cor- 
related ultracold atoms [FT;. Unfortunately, such pro- 
posal does not provide spatial resolution, and thus it 
cannot, e.g., discriminate between different antiferromag- 
netic quantum phases present in lattice models. Here we 
show that this important limitation can be overcome by 
using a standing wave probe laser configuration. Such 
a modification is crucial: it allows for spatially periodic 
QND coupling of light to atomic spins, and thus can re- 
veal spin correlations with the period of coupling. The 
quantum noise of the transmitted light then carries in- 
formation on the Fourier components of the spin density. 
Controlling the parameters of the optical lattice and of 
the probe light standing wave permits discrimination and 
characterization of various standard and "exotic" mag- 
netic orderings. 

In Fig. [I] we schematically depict various important 
lattice states, all of them having zero average of any 
component of the total spin. The paramagnetic state, 
Fig. [IJa), exhibits large (proportional to the total num- 
ber of atoms, iV a t) fluctuations of the total spin, whereas 
Fig.[IJb-d) show global singlets, i.e. states with total spin 
zero, which do not show any fluctuations. Obviously, 
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those most interesting strongly correlated global singlet 
states of Fig. [ljb-d) cannot be distinguished from each 
other by spatially homogeneous probing. Probing, how- 
ever, only every second atom in, e.g., the dimerized state 
Fig. |TJb) , i.e., measuring the total staggered magnetiza- 
tion of the state, will result in its correlations imprinted 
on the light. As we shall demonstrate below, such type of 
measurement offers the possibility to distinguish between 
various global singlets, with minimal disturbance to the 
system. 

Detection scheme The detection mechanism proposed 
here uses the off-resonant interaction of spin-F atoms 
(i.e., atoms having a 2F + 1 dimensional ground state 
manifold) with a polarized light beam propagating in z- 
direction (for a theoretical description of light-atom cou- 
pling cf. [HUE!])- The light pulse is characterized by the 
Stokes operators si, §2, S3 corresponding to the difference 
of the number of photons in the x and y linear polariza- 
tions, in the ±45° linear polarizations, and in the two 
circular polarizations, respectively. Let us first consider 
a one-dimensional (ID) sample of 2V a t atoms trapped in 
a ID optical lattice of period 7r/A:, oriented in z direction, 
with the n-th atom located at position z n = (n — l)ir/k, 
having the spin operator j(z„). The relevant dispersive 
part of the interaction Hamiltonian between light and 
atoms reads [THl HH] (see methods) 



H = ~KS 3 j 



(1) 



Here n is the coupling constant and J z cS is the z- 
componcnt of the effective collective atomic spin J = 
En c ™j( z ")- The coefficients c„ account for the modifi- 
cation of the atom-light coupling due to a spatial mod- 
ulation of the probe beam intensity (in a running wave 
configuration Cj = 1) and are the key parameters allowing 
for spatial resolution. As a particular case, we consider 
a standing wave configuration, as it is described in the 
caption of Fig. [2] Then, 

c n = c„(fc P , a) = 2 / dz cos 2 (k P (z - a))\w(z - z n )\ 2 , (2) 



with kp as the wavevector of the standing wave probe, 
and a its shift with respect to the optical lattice; w(z— z n ) 
is the usual Wannicr-type wave function of a single atom 
confined at z n in a deep optical lattice. 

We take the probe beam to be strongly polarized in 
the x direction so that Si (Si = f Sidt) fulfills (Si) = 
iVph/2 3> 1. This fact permits to introduce canonical 
quadrature operators X — S^/y^/Vph, P — Ssf-y/Np^, 
which satisfy [X, P] « i. Integrating the Heisenberg 
equation of motion for §2 shows that the effective collec- 
tive spin J° is imprinted on the X quadrature of light: 
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Since (Xj n ) = 0, the mean of this quadrature after pass- 
ing through the sample is directly proportional to the 



mean of the z-component of the effective atomic spin. 
Its variance contains the shot noise of the incoming probe 
pulse and the variance of J z cS . Assuming a coherent in- 
put beam; 
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The variance of the collective atomic spin can be ef- 
ficiently determined from the measurement on light if 
k = \/dr\ > 1 [16 . Here d = N at a/A denotes the res- 
onant column optical depth of the atomic sample and 
r\ = (NphcrT 2 ) / (A 2 A) is the probability of resonant exci- 
tation per atom by the probe, where a is the cross sec- 
tion on resonance for the probe transition, T the spon- 
taneous decay rate, A the detuning from resonance and 
A the cross section of the atomic ensemble illuminated 
by the probe. Since spontaneous emission destroys the 
spin state and the QND character of the coupling, the 
condition 77 <C 1 has to be fulfilled. To evaluate the de- 
coherence (and heating) caused by the probe, we need to 
estimate the effect of spontaneous emission on the vari- 
ance of the light quadrature. The induced probe decoher- 
ence is given by ((A J z ecoh ) 2 ) / FN at ~ r\ (in units of the 
coherent spin state noise). Therefore, the total measured 
noise reads 
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In order to measure the atomic spin fluctuations with 
the best possible accuracy^ we need to maximize the sig- 
nal to noise ratio [k, 2 ((A jf) 2 ) / FN at ]/[l/2 + k 2 tj}. For 
a coherent input beam, this contribution is maximized 
for 7] opt ss l/\/2d. Notice that 77 can be adjusted to rj opt , 
by choosing the appropriate detuning, intensity and du- 
ration of the laser probe. On the other hand, the optical 
depth of a sample of ultracold atoms, including e.g. a 
cubic lattice with 100 x 100 x 100 atoms can easily reach 
values of few hundreds. Thus, the value of q can be made 
much smaller than one |19j . Since the spin noise of inter- 
est (as shown later in Fig. 3) is of the order of unity, the 
additional noise provided by spontaneous emission does 
not, therefore, significantly modify the measured variance 
of the light quadrature. 

Furthermore, a squeezing of the X quadrature of the 
incoming probe before passing it through the sample such 
that ((AAj n ) 2 ) < 1/2, allows to improve the signal-to- 
noise ratio even further. Finally, let us mention that 
from the data recorded at the homodyne detector also 
higher order terms ((AA out ) m ) can be extracted. 

In order to achieve a spatial modulation of the atom 
light coupling, k p /k ^ 1 is necessary. The wavelength of 
the probe laser and thus kp is constrained by the near 
resonant condition for an efficient QND probing. Also 
the choice of the wavelength of the light forming the op- 
tical lattice is limited. Thus, tuning the wavelengths of 
the lasers only gives a restricted control over the ratio 
kp/k. An appropriate choice of the ratio can, however, 
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be achieved in various ways: (i) trivially by modifying k 
by changing the angle between the beams forming the op- 
tical lattice; (ii) more surprisingly, by shining the probe 
light at an angle 8 to the lattice. The effective Hamilto- 
nian then reads 

H = -ks 3 (cos 6J z cS - sin 6>J/ ff ) , (6) 

and the effective wavevector for a f D sample oriented in 
z-direction is kp = 2ttcos9/\p, with Xp the wavelength 
of the probe laser. Finally, (iii) a probe standing wave 
with variable kp can be obtained from crossing the probe 
laser beams at a variable angle. This last approach leads, 
however, to a different interacting Hamiltonian. 
Quantum antiferromagnets in ID Let us now illus- 
trate the power of our proposal, by applying it to concrete 
examples of strongly correlated states that can be real- 
ized with ultracold atoms. Particularly challenging are in 
this context various possible quantum antiferromagnetic 
states that lie in the center of interest of condensed mat- 
ter and even high energy physics [2T] |22] • We will analyze 
some of the states that appear as ground states (or ide- 
alizations thereof) of the generalized Heisenberg spin-1 
atomic chain [531 [Ml [55] ( see Methods). For some choice 
of parameters such systems arc in the, so called, Haldanc 
phase. This phase is gapped, and the ground state is 
well described by a matrix product state (MPS), that is 
both translationally and rotationally invariant. A promi- 
nent example of such states is the AKLT state [2U] (see 
Fig. [ljd)). In another regime of parameters the ground 
state exhibits the dimer, or Peierls, order and it repre- 
sents an example of the Valence Bond Solid (VBS) state. 
To a good approximation the ground state with an even 
number of sites is formed by placing a dimer on every odd 
bond (i.e., pairing the atoms on these bonds into singlets, 
see Fig. [ljb)). Such a state breaks the translational in- 
variance, but not the rotational symmetry. Yet for an- 
other choice of parameters, the system is in the critical 
phase, but exhibits some trimerization. The ideal t rimer 
state corresponds to concatenations of triples of neigh- 
bouring atoms forming a singlet, i.e. the state of zero 
total spin of the triple (Fig. [TJc)) . The result of apply- 
ing our detection method to the paramagnetic and these 
three types of global singlet states are shown in Fig. [3] 
While the presence of fluctuations at kp/k = sig- 
nals unambiguously the unpolarized paramagnetic state, 
the three distinct global singlets can be perfectly distin- 
guished cither by combining results for variable kp/k and 
a, or by fixing kp/k and varying solely a. For instance, 
at kp/k — 1/2 and a — 0, the added noise (in units of the 
shot noise of the probe light) for the dimerized, trimer- 
ized, and AKLT states is respectively 4k 2 /3, 8k 2 /9, and 
2k 2 . The above quoted values correspond the the ideal 
case in which the atoms are (^-localized on the minima 
of the optical lattice, while Fig. [3] shows the results for 
extended atomic Wannicr functions. This fact yields an 
overall decrease of the noise contribution towards larger 
kp/k ratios. Except for this decrease, the pattern is re- 
peated for larger values of kp/k. The fluctuations on the 



light quadrature using the atomic ground states of the 
spin-1 lattice systems obtained numerically [26], present 
a very close resemblance to the ones depicted in Fig. [3] 
Quantum antiferromagnets in 2D Our scheme is not 
limited to arrays of ID systems, ft can also detect and 
distinguish different antiferromagnetic ground states in 
lattice models in higher dimensions [21] [22]. As an il- 
lustration, we consider a Heisenberg spin-1/2 model on 
a square lattice with antiferromagnetic coupling. In this 
case it may happen that translational symmetry is broken 
and the ground state becomes the dimer VBS state (see 
Fig. |3Ja)), which exhibits long range dimer order. Also 
exotic spin liquid states which do not break any sym- 
metry and do not possess long range order can appear 
in an extended Heisenberg model. These states can be 
visualized as coherent superpositions of random dimer 
coverings, termed as Resonating Valence Bond (RVB) 
states (see Fig. Qb)). We consider a 2D lattice in the 
y — z plane, and the probe light strongly polarized along 
the ^-direction propagating in z direction. Assuming 
all atoms paired into dimers, we obtain that the addi- 
tional noise contribution to the outgoing X quadrature 
is proportional to the number of dimers oriented in the 
z-direction. A valence bond solid with all dimers ori- 
ented in the z-direction and a single valence bond state 
with randomly oriented dimers can then be easily distin- 
guished, since for kp/k = 1/2 and a = the additional 
noise is 2k 2 and k 2 , respectively. The antiferromagnetic 
Neel state appearing in the 2D Heisenberg model can be 
indeed detected with the method of Ref. [17] . 
Fermionic superfluidity Another application of our 
scheme is to detect superfluid phases in fermionic 
ultracold gases in a QND way. For the probing 
standing wave set in the z-direction the transmitted 
light carries information on the fluctuations of J z eff = 
J cPr cos 2 (kpz) J z (r). In a 2D or 3D system one can 
measure the fluctuations in J y and J x in the y and x 
direction, respectively, by changing the standing wave 
orientation. Similar quantities have been proposed and 
proved to detect superfluidity in ultracold Fermi gases 
[27] Hi US [301 EH]. Our scheme as described here can be 
specially interesting for a ID system as one can change 
the probing momentum kp with respect to the Fermi 
momentum by tilting the probing standing wave with re- 
spect to the system long axis as discussed above. Com- 
parison of these fluctuations for a paired and a non paired 
system reveals the existence of a pairing gap. Already in 
a homogeneous set-up one could measure the fluctuations 
in J z that, in a similar way as the structure factor, change 
dramatically across the BCS to BEC crossover [32] , 
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Methods 



coupling constant reads k = a\yj N^N^F /2. 

ID spin chains The strongly correlated states of spin- 1 
atoms on a ID lattice considered in this paper are ex- 
amples of ground states (or idealizations thereof) of the 
generalized spin-1 Hciscnberg Hamiltonian [23, 2H [2"5] 



Atom-light Hamiltonian The Hamiltonian of equa- 
tion is obtained from the off-resonant coupling to an 
atomic dipole transition. The excited atomic states can 
be adiabatically eliminated from the dipole-interaction 
Hamiltonian under the conditions described in the main 
text. The dispersive effects arising from the Stark shift 
of the atomic levels are described via an effective Hamil- 
tonian 



H = cos/3j(z„)-j(z„,)+sin/3 j(z n ) ■ j(z n ,) 

(n,n') 
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Here L is the length of the atomic sample and A its 
cross section overlapping with the probe light propagat- 
ing in the z-direction, and p = p{z) is the atomic density. 
The light is described via Stokes operators Sj = Si{z,t) 
{s± = §i ± 1S2) and 4> is the total photon density. The 
coefficients cij depend on the laser wavelength and the 
characteristics of the atomic transition. The first term, 
proportional to do, gives the (polarization independent) 
ac Stark shift. In the limit where the detuning is large 
compared to the hyperfine splitting of the excited state 
a2 ~ 0, and only the linear QND coupling between the 
Stokes operator and atomic spin remains. This is equiv- 
alently written in equation ([I]) as a coupling between the 
Stokes operator and the effective atomic spin component 
J z s . To proceed, we assume that J| ff is time indepen- 
dent. This is valid since: (i) it is conserved by the Hamil- 
tonian equation (JT|), and (ii) the measurement time for 
a pulsed probe is much shorter than the atomic spin dif- 
fusion time. The outgoing X quadrature of light is then 
given by Eq. (|| provided that a\F 2 NlJ2 < 1. The 



Due to the interplay between the bilinear and the bi- 
quadratic interaction (parametrized by 0), this Hamilto- 
nian presents three antiferromagnctic quantum phases: 
dimerized /3 € (— 3tt/4, — 7r/4), Haldane f3 E (— 7r/4, 7r/4), 
and critical (3 € [n/4,n/2). The representative point 
of the dimerized phase is found at j3 = —ir/2 where 
'the two-site ground state is a singlet (dimer, [531 l2~4"]l: 
|s)i,2 = (|+1>i|-1)2+|-1)i|+1>2-|0)i|0) 2 )/V3. Here we 
denote by \m z = ±l,0) n the eigenstates of j z (z n ). This 
ground state is analytical only in the thermodynamic 
limit or for periodic boundary conditions. In a finite 
chain with an even number of sites, however, the ground 
state is close to a concatenation of spin-1 singlets (see 
Fig.[IJb)). The representative state of the Haldane phase 
is found at /3 = arctan(l/3), where the exact ground 
state corresponds to the AKLT state (see Fig.[ljd)). The 
critical phase is less well understood, though known to 
have period three correlations. At j3 — ir/4 the three- 
site ground state is a singlet (trimer), and a caricature 
of the ground state for a larger number of sites is given 
by a concatenation of trimcrs (cf. Fig.Qc)). These three 
isotropic ground states are eigenstates of any component 
of the total spin. Thus for kp = k where all atoms couple 
equally to the light, no additional noise is imprinted on 
the X quadrature. For kp ^ k, in general these states 
are not eigenstates of J z cS (although (J z cti ) = 0) and as 
we demonstrate in Fig. [ljc), the fluctuations permit to 
discriminate them unambiguously. 
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FIG. 1: Antiferromagnetic states of spin-1 lattice sys- 
tems, (a) 3D cubic lattice with a paramagnetic state of un- 
polarized atoms; (b) dimerized state with pairs of neighbour- 
ing atoms forming singlets; (c) trimerized state with triples 
of neighbouring atoms forming singlets; (d) AKLT (Affleck- 
Kennedy-Lieb-Tasaky) state obtained from a concatenation 
of spin- 1/2 singlets by projecting pairs of spins from different 
bonds into the subspace of total spin-1 20 . 




FIG. 2: Schematic experimental setup. A strong laser 
beam initially polarized in x-direction is impinging on a 99/1 
beamsplitter. The transmitted attenuated part of this probe 
is propagating through the sample and reflected off a mirror, 
such that a standing wave with wavevector kp in z direction 
is formed at the position of the atomic sample. After the 
second pass, the laser beam is outcoupled to a homodyne 
detector, where (S-i) is recorded. The atoms are trapped in 
an optical lattice with wavevector k. Its relative displacement 
with respect to the probe standing wave can be changed by 
either moving the end mirror, or by introducing a phase shift 
between the counterpropagating laser beams. 
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FIG. 3: Detection of antiferromagnetic states of spin-1 lattice systems. Fluctuations (e = ((AX out ) 2 } — 1/2) imprinted 
on the X quadrature of the probe light beam transmitted through a sample of spin-1 atoms in a one-dimensional lattice for (a) 
unpolarized paramagnetic, (b) dimerized, (c) trimerized, and (d) AKLT states. While the presence of fluctuations at kp/k = 
signals the unpolarized paramagnetic state, the other three global singlets phases can be unambiguously distinguished at 
different values of the ratio kp/k and/or the shift a (see text for details). 



FIG. 4: Two-dimensional states of dimers. (a) Valence 
Bond Solid arranged in a regular crystal, (b) Valence Bond 
Solid consisting in a random covering of dimers. 



